Preparing the bound instance of quantum entanglement 
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Among the possibly most intriguing aspects of quantum entanglement is that it comes in "free" and "bound" 
instances. Bound entangled states require entangled states in preparation but, once realized, no free entangle- 
ment and therefore no pure maximally entangled pairs can locally be regained. Their existence hence certifies 
an intrinsic irreversibility of entanglement in nature and suggests a connection with thermodynamics. In this 
work, we present a first experimental unconditional preparation and detection of a bound entangled state of 
light. We consider continuous-variable entanglement, use convex optimization to identify regimes rendering 
its bound character well certifiable, and realize an experiment that continuously produced a distributed bound 
entangled state with an extraordinary and unprecedented significance of more than ten standard deviations away 
from both separability and distillability. Our results show that the approach chosen allows for the efficient and 
precise preparation of multi-mode entangled states of light with various applications in quantum information, 
quantum state engineering and high precision metrology. 



PACS numbers: 03.67.Mn, 03.67.-a. 03.65.Ta, 03.65.Ud 

The preparation of complex multi-mode entangled states of 
light distributed to two or more parties is a necessary starting 
point for applications in quantum information processing 1 1 1 
[21 as well as for fundamental physics research. An aggres- 
sively pursued example of the latter is the preparation of the 
bound instance of entanglement, a type of entanglement that 
can only exist in higher-dimensional or multi-mode quantum 
states |3|. Bound entanglement is fundamentally interesting 
since, in contrast to "free" entanglement, it can not be dis- 
tilled to form fewer copies of more strongly entangled pure 
states [3 1 by any local device allowed by the rules of quantum 
mechanics. This irreversible character has triggered entire 
theoretical research programmes |4|, in particular by linking 
entanglement theory to a thermodynamical picture, with this 
irreversibility reminiscent of — ^but being inequivalent with — 
the second law of thermodynamics [5J. In order to investigate 
such connections both new theoretical as well as experimental 
means of constructing multi-mode states must be innovated. 

In recent years, great progress in information process- 
ing, metrology and fundamental research has actually been 
achieved in the photon counting (discrete variable, DV) 
regime using postselection |1|. States of light are the opti- 
mal systems for entanglement distribution because they prop- 
agate fast and can preserve their coherence over long dis- 
tances. Postselection means that the measurement outcome 
of the detectors which characterizes the quantum state is also 
used to select the state, conditioned on certain measurement 
outcomes. In such an approach, conditional applications are 
possible, however, an unconditional application of the states 
in downstream experiments is conceptually not possible. An- 
other limitation that any postselected architecture will even- 
tually face is that without challenging prescriptions of mea- 
surement, quantum memories and conditional feedforward, 
the preparation (post-selection) efficiency will exponentially 



decay with an increasing number of modes. In parallel to 
postselected architectures of light, unconditional platforms for 
research in quantum information have been developed which 
build on the detection of position and momentum like vari- 
ables having a continuous spectrum and a Gaussian statistics. 
In such platforms the preparation efficiency of one mode is 
identical to the preparation efficiency of N modes. In the 
past, this continuous variable (CV) platform has been used to 
demonstrate the Einstein-Podolsky-Rosen (EPR) paradox Q 
and unconditional quantum teleportation |7|. Recently, the 
CV platform has been extended to investigate multimode en- 
tangled states 1 8 1; however, the significance of their nonclas- 
sical properties have typically been smaller compared to their 
postselected counterparts. 

In this work, we demonstrate the continuous unconditional 
preparation of one of the rarest types of multi-mode entan- 
gled states - bipartite bound entangled states - using the CV 
platform. The property of bound entanglement is verified by 
four downstream balanced homodyne detectors with a detec- 
tion efficiency of almost unity. Alternatively, our setup can 
make available bound entangled states for any downstream 
application. The bound entanglement is generated with un- 
precedented significance, i.e., with state preparation error bars 
small with respect to the distance to the free entanglement 
regime and with respect to the distance to the separability 
regime. Our result is achieved by the convex optimization 
of state preparation parameters, and by introducing the exper- 
imental techniques of single-sideband quantum state control 
and classical generation of hot squeezed states. 

The first ever generation of bound entangled states was 
claimed in 2009 |9|. This work used photon counting and 
postselection, however, the data presented did not support this 
claim, an issue which has been addressed in a comment, see 
Ref. iflOl . In Ref. ifTTl a DV nuclear magnetic resonance state 
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whose density matrix has a small contribution of bound en- 
tanglement has been observed. Such a state has been called 
a "pseudo-bound entangled state". Very recently, the actual 
first bound entangled states have been generated in two ex- 
periments, both on the basis of discrete variables. In Ref |i,12J 
bipartite bound entangled states of trapped ions have been ver- 
ified by the unconditional detection of resonance fluorescence. 
In Ref. 1 1 3 1 the first bound entangled states of light have been 
generated, albeit of multipartite and not of bipartite nature. 
Similar to Ref |9|, photon counting and postselection have 
been used. An unconditional application of the distributed en- 
tanglement in a downstream experiment is hence not possible. 
This is now made possible in our work, with a significance of 
bound entanglement that has not been achieved using postse- 
lection. 

Our theoretical search for CV Gaussian bound entangled 
states of light begins with three (non-pure) squeezed input 
modes and a vacuum mode overlapped on four beam split- 
ters acting as phase-gates. This yields several independent 
parameters to be chosen that includes three pairs of quadra- 
ture variances and the splitting ratios and the relative phases 
of the phase-gates. Additional vacuum contributions due to 
optical losses at different locations in the experiment have to 
be considered as well. As it turns out, bound entanglement 
is extremely rare in this multi-dimensional parameter space. 
Hence, to theoretically identify suitable regimes for experi- 
mental certification is a challenging task: Known examples of 
CV bound entangled states, including those of Ref 1 14|, will 
have both free entangled and separable states very nearby. Op- 
timal entanglement witnesses can be efficiently constructed 
for Gaussian states [16], yet to maximize the distance of an 
optimal hyperplane separating separable states to the bound- 
ary of non-distillable states — hence maximizing robustness of 
a preparation — is a non-convex difficult problem. What is 
more, a reasonable compromise with the preparation com- 
plexity has to be found, with a surprisingly simple feasible 
scheme being shown in Fig.[T] 

We now present the measures required for verifying the 
presence of bound entanglement. Since the studied states 
are Gaussian they are fully described by their first — which 
will not play a role here — and second moments, specified 
by the covariance matrix of a state p lfT7HT9l . We define 
a set of quadratures for each optical mode given by Xj ~ 
{cLj + aj)/2^/-^ and pj = —i{aj — aj)/2^/^ where aj,a| are 
the annihilation and creation operators, respectively. Collect- 
ing these 2n coordinates in a vector O = {xi ,pi, . . . , £„, p„), 
we can write the commutation relations as [Oj, Ok] = io'j.k, 
where h — 1 and is a matrix a often known as symplectic 
matrix. The second moments are embodied in the 2n x 2n 
covariance matrix 

7j- = 2Re tr (/}(6j - ) (O^ - 4 )) , ( 1 ) 

with dj = tr(pOj), giving rise to a real-valued symmetric 
matrix 7, see supplementary material. 

Verification of bipartite bound entanglement requires show- 



ing that the state is entangled (inseparable) with respect to a 
bipartition of the modes and that the state remains positive un- 
der partial transposition [|3][l4l[T5l proving that the state is not 
distillable. The state is said to be entangled if no physical co- 
variance matrices •ja and jb exist of states in modes A and B, 
respectively, so real matrices satisfying 7a, 7s > —if, such 
that IIT4I 7 > 7A ® 7s- This idea suggests a natural entangle- 
ment measure ll20l for Gaussian states, defined as the solution 
of 

E{j) = 1 - max X (2) 
7 > 7A ffi 7s, 7A, 7s > -ixa. 

E{'-f) > indeed implies that the state is entangled. The 
above problem is known as a semi-definite program, a convex 
optimization problem that can efficiently be solved. 

Non-distillability can be tested by evaluating the par- 
tial transposition of a state lO) which physically reflects 
time reversal. For covariance matrices, partial transposi- 
tion amounts to changing the sign of momentum coordinates 
or by applying the operation 7^^ = M^M, where M — 
(1, 1, 1, 1, 1, —1, 1, —1), with a —1 in all momentum coordi- 
nates belonging to B. A covariance matrix 7 is said to be PPT 
if its partial transpose is again a legitimate covariance matrix, 
or equivalently, 7^" + icr > 0. A measure as to the quanti- 
tative extent a state is PPT can be taken to be the minimum 
eigenvalue of this matrix, 

P(7) = mineig(7^ + icr). (3) 

The continuity of the eigenvalues with respect to variations in 
the matrix are enough to guarantee that the measure is mean- 
ingful. A strictly positive value of P{'^) unambiguously certi- 
fies that the state is not distillable. 

Based on our theoretical parameter search our final exper- 
imental setup is realized as shown in Fig. [T] In total three 
optical parameter amplifiers (OPAs), three phase-gates, con- 
sisting of a beam splitter and a piezo mounted mirror, and a 
vacuum mode are utilized as the base setup. The four homo- 
dyne detectors are only necessary for the verification of bound 
entanglement but not for its preparation. We set our OPAs to 
produce the minimum and maximum vacuum noise normal- 
ized vai-iances to be: (2.0, 3.46) from OPAi, (0.54, 5.16) from 
OPA2 and finaUy (0.63, 2.54) from OPA3. The phase-gates 
were set to 0i = 90°, 02 = 41° and (f)^ = 140°, respectively. 
For further details, see Ref. ETI . 

The first OPA produces a classically squeezed (thermal) 
state we refer to as hot squeezing. It manifests a non-uniform 
stationary noise distribution amongst its two quadratures with- 
out having the smallest quadrature fall below the vacuum 
noise level. Hot squeezing is generated when, for example, 
two amplitude squeezed modes of different squeezing factors 
are overlapped on a 50/50 beam splitter with a relative phase 
of 90°, thereby producing a two-mode squeezed state, but then 
one of the output modes is discarded. Without the presence 
of hot squeezing, we could not find any possibility to pro- 
duce continuous-variable bound entanglement (between four 
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FIG. 1. Experimental setup: The experiment is composed of three 
optical parametric amplifiers (OPAi-a), three actively controlled 
piezo mounted mirrors forming phase-gates (PGl — 3) and four ho- 
modyne detectors which are independent of the preparation. The in- 
set shows the construction of an OPA as a non-linear crystal inside a 
resonator producing a spatial TEMoo (transverse electro-magnetic) 
mode. The bound entangled state is obtained through the bipartite 
splitting such that Alice and Bob each possess two of the four modes. 



modes), using the methods described in Ref. |2T|. Indeed, the 
passive optics following the sources can no longer alter the 
eigenvalues of 0-7, which also define the degrees of squeezing 
and mixedness of the state. Hot squeezing therefore appears 
to give rise to a necessary ingredient of quantum and classical 
correlations in order to create robust bound entangled states. 
We demonstrate that the same state can also be prepared in a 
purely classical way by applying a local random displacement 
on the phase quadrature of a vacuum mode while parametri- 
cally amplifying the state's amplitude quadrature. The station- 
ary random phase modulation is produced by using an EOM 
(electro-optical modulator) driven with the output from a ho- 
modyne detector measuring shot noise. The amplitude mod- 
ulation is generated by operating OPAi in Fig. [T] in amplifi- 
cation mode, effectively anti-squeezing the amplitude quadra- 
ture and deamplifying the thermal noise phase quadrature of 
the input state. In principle the random amplitude noise of 
the first input mode can also be provided by a second homo- 
dyne detector and an amplitude modulator, thereby replacing 
the parametric OPAi device. We note that pseudo-random 
numbers could be insufficient in this scheme since they could 
introduce artificial correlations and a non-stationary noise into 
the final state. 

In order to hit the tiny regions in parameter space where 
bound entanglement does exist we introduce to our setup a 
new technique for precisely controlling phase-gates at arbi- 
trary angles. This method relies on an optical single-sideband 
scheme (see supplementary material) that can be used to ar- 
bitrarily and independently set the working point of both a 
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FIG. 2. Experimental results: The state measured after 4 million sets 
of raw quadrature data points yields the entanglement E and non- 
distillability P indicated by the red cross. Other lO'' points are ob- 
tained by bootstrapping the original 4 million data points and show 
that we are 16a away from separability and 46(t away from distill- 
ability. In the inset we depict the minimum eigenvalue of 7 + icr of 
each of the 10* bootstrapped correlation matrices, showing that they 
are significantly far away from the boundary of covariance matrices 
allowed by the uncertainty principle. 



phase-gate network and multiple homodyne detectors. This 
scheme reduces setting the relative phase between interfering 
modes to selecting the electronic demodulation phase used in 
the control loop. A portion of the light leaving the phase- 
gates, PGl-3 in Fig. [T] is redirected to control photodetec- 
tors. We are able to derive a strong error-signal by tapping 
only 1 /xW of power corresponding to no more than 1 % of the 
signal mode's optical power. For applications where delicate 
quantum states must remain free from losses our method pro- 
vides a means by which they can still be used for controlled 
interference without significant vacuum contribution due to 
loss. 

The four balanced homodyne detectors are used for the full 
tomographic reconstruction of the covariance matrix. The 
results of the reconstruction are used to evaluate two char- 
acteristics of the state; namely, its entanglement E (2) and 
its PPTness P Eq. ([3]). In order to build the statistics of 
these characteristics we first continuously recorded 4 million 
data points from the amplitude and phase quadratures of each 
mode. Using the bootstrapping method, we then randomly 
sampled from the total 4 million points, with uniform distri- 
bution, points that were different, and produced a series of 
covariance matrices from which the entanglement, PPT and 
physical properties were calculated. Our results are repre- 
sented in Fig. |2]by the black points. The cross corresponds 
to the average state inferred from the total data set. The ab- 
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scissa of Fig.|2]is the PPTness and the ordinate the entangle- 
ment. By projecting the scatter plot onto the respective axes 
we calculate a significance of 46(7 away from being distill- 
able, i.e., P(7) < and 16(t away from being separable, i.e., 

(7) < 0. To demonstrate that the generated state is not close 
to the boundary of state space (and to confirm its physicality), 
the smallest eigenvalue of 7 + jct is also depicted: This is 
shown in the inset as a histogram. The fact that it is more 
than 50(7 away from being unphysical can be seen as an in- 
dication of the fact that our setup was stable over the entire 
measurement time and that our measured data exhibited little 
statistical uncertainty. 

The fact that the involved states are Gaussian to a very high 
statistical significance is carefully checked by statistical meth- 
ods (see supplementary material |21|): Indeed, we have not 
only estimated the second moments from the time series data 
from homodyning measurements, but in fact all moments. 
Based on these data, we have computed Q-Q-plots for mea- 
sured distributions against perfectly Gaussian ones, showing 
a remarkable coincidence and confirming the precisely Gaus- 
sian character of the state. What is more, from a bound of the 
mean energy of the state and of its Gaussian character, one 
can derive rigorous bounds to the distillable entanglement of 
the state, confirming an at most negligible amount of distill- 
able entanglement. Again, for details, see the supplementary 
material [21 j . 

Our results present the first unconditional preparation of 
bound entangled states of a physical system characterized by 
(continuous) position/momentum-like variables. With respect 
to systems composed of light, we demonstrate the first un- 
conditional preparation of bound entanglement, and achieve 
an unprecedented significance of its features. Independent of 
any postselection, our platforms allows for the distribution of 
the entangled states. As other states of light our bound entan- 
gled states can be distributed to remote parties, which might 
be kilometers apart using optical fibers ll22l . The decoher- 
ence on bound entangled states due to photon loss and phase 
noise [23] and the ineffectiveness of distillation schemes |24| 
can be tested, as well as the applicability of thermodynamic al 
pictures of entanglement be studied experimentally. 

Our results clearly exemplify the potential of the continu- 
ous variable platform for the precise engineering of complex 
multi-mode states of light. We underline that using this plat- 
form the state preparation efficiency does not depend on the 
number of entangled modes. That is to say, detecting, for ex- 
ample, one squeezed mode with one homodyne detector has 
exactly the same efficiency as detecting squeezed states 
with A^ homodyne detectors simultaneously. Furthermore, we 
estimate our total quantum detection efficient to be between 
90-95% being already considered in the preparation of bound 
entanglement. Alternatively, this loss could be mapped di- 
rectly onto the measured state by inclusion of neutral density 
filters, and verification with perfect detectors would reveal the 
same statistics as depicted in Fig.|2] 

We believe that the precise and unconditional preparation 
of (bi-partite) bound entangled states of light demonstrated 



uplifts the theoretical and experimental research on the link 
between entanglement theory and statistical physics. From 
a more general and also technological perspective, the high 
efficiency and the high degree of control in multimode 
quantum state preparation achieved certainly promotes the 
application of the unconditional continuous variable platform 
for the preparation of quantum states of light for fundamental 
research as well as quantum metrology. 

This work has been supported by the EU (QESSENCE, 
MINOS, COMPAS), the EURYI, the grant UNAM-PAPIIT 
IN 11 73 10 and by the Centre for Quantum Engineering and 
Space-Time Research, QUEST. We acknowledge discussions 
with P. Hyllus and M. Lewenstein at an early stage of this 
project. 



Appendix A: Heisenberg uncertainty and entanglement criteria 

Explicitly, for n modes the symplectic matrix a reads as 

3=1 ^ ' 

The Heisenberg uncertainty relation, expressed in terms of 
the covariance matrix [ 17], is given by 

7 + icr > 0. (5) 

Such operator valued inequalities A > B for Hermitian A 
and B always refer to operator ordering, meaning that the real 
eigenvalues of A — B are non-negative. The entanglement 
measure E for covariance matrices defined in the main text in- 
deed indicates entanglement in states ll20l . and for two modes 
this is essentially nothing but the familiar negativity ||25]| . 

In the discussion of the main text we show that the spectrum 
of 7 + icr is bounded from below by e > 0, hence manifesting 
the Heisenberg uncertainty principle. It is worth mentioning 
that this also means that the smallest symplectic eigenvalue 
51(7) of 7 is bounded away from 1. In the experiment, we 
also test whether the reconstructed covariance matrix satisfies 
the Heisenberg uncertainty relation as this is a test if the ma- 
trix corresponds to a physical state. Unphysical states might 
occur if the error bars of the quantum state preparation or the 
tomographic characterization are too large. 



Identifying robust bound entangled states 

The relative volume of bound entangled states compared to 
all states is very small under every reasonable measure, and 
any verification as pursued here necessarily requires a careful 
analysis as to what parameter set is most suitable. In this sub- 
section, we report techniques that have been used to identify 
regimes of robust bound entangled states. 
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Consider a general 4 mode correlation matrix, expressed in 
bipartite normal form ETl : 
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The 16 parameters can be seen as describing a manifold in 
R^^. We sample uniformly the hypercube [—1/2, 1/2]^^^ un- 
til we get a bound entangled state. The evaluation of the P 
measure for a given covariance matrix 7 amounts to solving 
an eigenvalue problem, that of the degree of entanglement E 
to solving a semi-definite problem. In practice, the latter prob- 
lem can also be performed in the dual space of witnesses that 
are quadratic polynomials in the quadratures, as explored in 
Ref |16|. Once an instance is found we construct a random 
walk in order to improve the robustness, by (i) displacements 
in the direction of the axis by a small amount 5X ~ 0.01 and 
(ii) rotations, by the same angle, in each of the 16 x 15/2 
two dimensional Cartesian planes; the new covariance matrix 
is accepted if the new corresponding state has a larger degree 
of entanglement and is more significantly a PPT state as mea- 
sured by E and P, respectively. The most suitable state found 
(after several hundred hours of computer time), as quantified 
by the biggest value of min{i?(7), P(7)}, is characterized by 
an entanglement value of E{^) = 0.054 and P(7) = 0.132, 
giving an idea of the limiting values that one can achieve. 

However, experimentally it is too expensive to engineer a 
state with an arbitrary correlation matrix. We thus construct a 
circuit which, starting from a product of noisy Gaussian single 
mode states, can produce bound entangled states, but is simple 
enough to be producible in the lab with available technology. 
A (non-unique) example of such a circuit is plotted in Fig. 1. 
The resulting scheme is a result of a variation within the given 
parameterized family of circuits - again using a random walk 
approach as described above, but now on the physically fea- 
sible set of covariance matrices by parametrizing each of the 
optical components - maximizing the statistical significance 
of being bound entangled by running semi-definite problems 
in each step. Afterwards we filter the results allowing only 
those which require achievable values of squeezing at the in- 
put and which only require a single mode with hot squeezing, 
as this is also a precious resource that, at the moment, can 
only be input in a single mode. Within the resulting states 
we choose the most robust according to the aforementioned 
criteria. 



Appendix B: Details of the experiment 

The three OPAs used to produce the underlying quadrature 
squeezing at sideband frequency of 6.4 MHz were constructed 



from a type I non-critically phase-matched MgO:LiNb03 
crystal inside a standing wave resonator, similar to the de- 
sign that previously has been used in Ref. Il26l . They were 
pumped with approximately 100 mW of green light at 532 nm 
each resulting in a classical gain of about 5. The length of 
the OPA cavity as well as the phase of the second harmonic 
pump beam were controlled using radio-frequency modula- 
tion/demodulation techniques. 

Balanced homodyne detection was performed on each of 
the four modes in order to reconstruct the 8x8 covariance ma- 
trix. The optical local oscillator was filtered through a three 
mirror ring cavity operated in high finesse mode resulting in 
a linewidth of 55 kHz. The detector difference currents were 
electronically mixed with a 6.4 MHz local oscillator and low- 
pass filtered with a 400 kHz bandwidth. The dark noise sepa- 
ration from shot noise was measured to be more than 10 dB for 
each detector. The raw data was acquired using a 14 bit Na- 
tional Instruments DAQ-card and in total eight measurement 
settings including the shot noise measurement were required 
in order to reconstruct the covariance matrix. 

The hot squeezed states were generated by randomly phase 
modulating the control beam used to set the length of the 
OPA cavity at the squeezing sideband frequency, 6.4 MHz, 
and locking the OPA cavity in amplification. This produces 
phase squeezed states whose smallest quadrature can be con- 
trolled by varying the strength of the random noise modulated 
on the control field and whose amplitude quadrature is con- 
trolled by the degree of classical gain. 

The single-sideband was generated by overlapping the out- 
put of a second laser operating at around 1064 nm with the 
bright output of OPAl. The beams were phase-locked at a 
beat frequency of 15 MHz resulting in a field that corresponds 
to both a phase and amplitude modulation. The beat was de- 
tected by directing approximately 1% of the phase-gate out- 
puts to photodetectors placed behind the phase-gates as well 
as in each homodyne detector. The relative phase between the 
carriers at both the phase-gates and the homodyne detectors 
could then be set to an arbitrary phase simply by changing 
the demodulation phase of the electronic local oscillator We 
estimate a phase sensitivity at each phase-gate to be approxi- 
mately 2deg. 



Appendix C: Discussion of the Gaussian character of the state 

The findings of our work show with a remarkable and un- 
precedented statistical significance that the second moments 
are certified to be those of Gaussian bound entangled states. 
What is more, a perfectly Gaussian statistics is clearly ex- 
pected due to the underlying physical mechanisms of our 
setup, i.e., parametric amplification of vacuum states, their su- 
perposition on beam splitters, time-independent linear losses 
and phase shifts, and finally balanced homodyne detection. In 
this subsection, however, we highlight the extent to which the 
states are indeed Gaussian states, and carefully discuss some 
theoretical issues associated with the certification of bound 
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entanglement for infinite-dimensional quantum systems. 

In the course of our work we have taken time series data 
from the homodyning measurements and have not only es- 
timated the second moments (leading to the covariance ma- 
trix specified in the main text), but in fact all moments. For 
each moment, we have tested for the Gaussian character of 
the state. Indeed, such studies of Gaussianity are interesting 
in their own right. More precisely, for our experiment, we 
have computed Q-Q-plots for a measured distribution against 
a Gaussian one, comparing these two probability distributions 
by plotting their quantiles against each other. Fig. [3] shows 
a representative sample Q-Q-plot of this kind. In each case 
the Gaussian character is verified to a remarkable level of sta- 
tistical significance. We also performed a test which also 
confirmed normality. Hence, one can conclude that the pre- 
pared is indeed Gaussian to a very high degree of accuracy. 
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Theoretical Quantiles 

FIG. 3. Q-Q-plot of the sample quantiles of a measured 4 million 
point data set versus the theoretical one for a Gaussian distribution 
based on the measured second moments (red dots). For comparison 
we also give a Q-Q-plot by plotting quasi-random numbers from a 
perfect Gaussian distribution (black dots). Such a Q-Q-plot depicts 
the q-th quantile of one distribution against the g-th quantile of the 
other. The diagrams show that the prepared states of our work are 
Gaussian to a very high level of significance. 

Having said that, strictly speaking, one may argue that in 
infinite-dimensional Hilbert spaces, the free-entangled states 
are (trace-norm) dense in state space |28j,,29J. Hence, in the 
vicinity of every bound entangled Gaussian state one can find 
a free entangled state that is operationally indistinguishable. 
This is as such no surprise at all: In the same way, it is true 
that arbitrarily close to any separable state there is a free en- 
tangled state for continuous-variable systems. However, in 
quantitative terms, the degree of free entanglement will indeed 



be negligible ll30l . That is to say, the measured mean energy 
and the closeness to Gaussian states readily give rise to rig- 
orous bounds to the distillable entanglement of the unknown 
state. Therefore, even in quantitative terms, one can falsify the 
proposition that a significant distillable entanglement is found 
- again strengthening the claim made in the main text. 
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